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INDUCED PRESSURE AND TEMPERATURE GRADIENTS IN A 
PLASMA AS A RESULT OF WEAK ELECTRIC FIELDS 

By Willard E. Meador 
Langley Research Center 

SUMMARY 

Computations are presented of the induced pressure and temperature gradients 
resulting from the application of weak electric fields to plasmas. Two methods are 
employed to solve an approximate model of the corresponding electron Boltzmann equa- 
tion: (1) the direct approach consisting of an exact solution through quadratic terms in 
the electric field, and (2) Everett's technique for closing out the macroscopic equations 
of change with Grad's 13-moment velocity distribution function. Comparisons between 
the two sets of results indicate substantial differences for most interparticle interaction 
potentials, especially with regard to the induced temperature gradient; hence, the Grad 
13-moment approximation may not accurately describe some of the fundamental aspects 
of plasma transport phenomena. 


INTRODUCTION 

A general assumption of many authors (see ref. 1 for listing) over the past 12 years 
is that Grad's 13 -moment velocity distribution function accurately describes nonequilib- 
rium plasmas when it is used to close out the macroscopic equations of change. Since, 
however, the Grad approximation is specifically an expansion (ref. 1) in the eigenfunctions 
of Boltzmann's binary elastic collision operator for Maxwellian molecules (that is, inverse 
fourth-power interparticle interaction potentials), a rapid convergence rate is not guaran- 
teed for more general potentials; consequently, additional justification is required for 
most plasmas of practical interest. Previous studies (refs. 2 and 3) of the convergence 
properties of Grad-like calculations have shown that although the residual heat flux, the 
entropy density, the collisional production rate of entropy, and second-order contributions 
to the electron pressure tensor are predicted with reasonable accuracy (within 26 percent) 
for spatially homogeneous systems, third-order contributions to the electron current den- 
sity and the phase angles by which the current density and the heat flux lag an applied 
oscillating electric field are off by factors up to 4 in the limits of very hard or very soft 
particles. 



The purpose of the present research is to extend the aforementioned convergence 
studies to spatially inhomogeneous plasmas. In particular, the radial pressure and tem- 
perature gradients induced in a plasma by the application of an axial electric field are 
computed by using two methods to solve an approximate model of the corresponding elec- 
tron Boltzmann equation. These methods are (1) the direct approach consisting of an 
exact solution through quadratic terms in the electric field, and (2) Everett's technique 
(ref. 4) for closing out the macroscopic equations of change with Grad's 13-moment veloc- 
ity distribution function. A wide range of effective interparticle interaction potentials is 
employed in both sets of calculations with the following results: Whereas significant dif- 
ferences in the pressure gradients occur only for particles which are much harder or 
much softer than the Maxwellian variety, the Everett-like predictions of the temperature 
profiles are very poor even for Maxwellian molecules. Hence, a feature which was never 
observed in the previous studies appears in the present investigation. 

SYMBOLS 


A ratio of collision integrals 

b impact parameter 

c e electron particle velocity relative to laboratory frame of reference 

Cj,C 2 ,C 3 transport coefficients 

e magnitude of electron charge 

E applied electric field 

f e electron velocity distribution function 

f^°) Maxwellian contribution to electron distribution function 

i,j indices 

i,j,k unit vector in x-, y-, and z-direction, respectively, of Cartesian coordinate 

system 

j electron current density 


2 



k 


m 


e 


n e 


n h 


Pe 



q 

Q 


r 


t 


T e 


T* 

A e 


u 


Boltzmann’s constant 
electron particle mass 
electron number density 

spatially homogeneous contribution to electron number density 
number density of heavy particles 

electron partial pressure relative to laboratory frame of reference 
spatially homogeneous contribution to electron partial pressure 
electron partial pressure relative to electron frame of reference 

traceless electron pressure tensor relative to laboratory frame of reference 

traceless electron pressure tensor relative to electron frame of reference 
heat flux 

dimensionless residual heat flux 
radial distance 

integral defined by equation (12) 
time 

electron temperature relative to laboratory frame of reference 
spatially homogeneous contribution to electron temperature 
electron temperature relative to electron frame of reference 
reduced electron particle velocity relative to electron frame of reference 
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U unit tensor 

v e electron diffusion velocity 

x integration variable 

x,y,z Cartesian coordinates 

X dimensionless pressure derivative 

X e electron body force per unit mass 

Y dimensionless temperature derivative 

/3q reduced electron diffusion velocity 

/3j reduced heat flux relative to laboratory frame of reference 

reduced heat flux relative to electron frame of reference 
jSj general nondimensional flux 

y reduced electron particle velocity relative to laboratory frame of reference 

e azimuthal angle for collisions 

£ effective interparticle interaction parameter 

a electrical conductivity 

t collision time 

<pi ith-order electron perturbation function 

X scattering or deflection angle 

Notation: 

Primed quantities in collision integrals denote conditions after a collision, as opposed to 
unprimed quantities which denote conditions before a collision. 
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(8f e /8t) c 



collisional time derivative of f e 
average over velocity space 

METHODS OF ANALYSIS 


The approximate model of the electron Boltzmann equation employed in the present 
research for both the exact and the Everett-type solutions is one in which Meador's sim- 
plified collision model (ref. 5) replaces the more precise collision integrals of Chapman 
and Enskog (ref. 6). The basic assumption of the collision model is that an effective 
interparticle interaction parameter | , which appears formally as the exponent in a sin- 
gle type of inverse-power electron — heavy-particle interaction potential, can be chosen 
semiempirically to include the effects of electron-electron encounters and multiple heavy 
species. The advantage of the model is the ease with which the corresponding integrodif- 
ferential kinetic equation can be reduced to differential form and then solved. 

Several examples of the reliability of this approach are given in reference 5 for 
plasmas involving electron-electron collisions; in addition, the parameter £ has a pre- 
cise physical meaning for Lorentz plasmas — that is, slightly ionized gases or full ioniza- 
tion with large ionic charges. Consequently, the present results and comparisons should 
provide strong implications about the general accuracy of the Grad 13 -moment 
approximation. 

The first step in the exact method of solution is the assumption that the applied 
electric field E is sufficiently weak to permit a power series expansion of the electron 
distribution function in the form 



where fg°^ is the Maxwellian contribution and </q is proportional to the ith power of 
E. This expression is substituted into the approximate Boltzmann equation and the coef- 
ficients of like powers of E are equated to obtain a sequence of differential relations, 
each of which is solved exactly for the corresponding dq. If solutions are obtained in 
this manner through $ n , the resulting electron distribution function is called nth order 
and is said to be exact through terms containing E n . Only second-order functions are 
considered in the present research. 

An important feature of the exact method is that some members of the aforemen- 
tioned sequence of differential relations will contain terms which do not appear in the 
analogous perturbation equations of Chapman and Enskog (ref. 6); in particular, may 
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occur in the ith-order equation in terms other than those derived from the collision inte- 
grals. Some precedence for this character is found in reference 6 in connection with 
magnetic forces, the inference now being that more general causes may exist. In any 
case, the present equations never contain fewer contributions than those of Chapman and 
Enskog and thus represent a generalization of that technique. 

Unlike the exact nth-order method, which gives the complete contribution to f e 
through E n terms, the Everett technique only partially describes the nth-order distri- 
bution function because the value of each term is affected by the number of moments 
employed. An unquestioned application of the Grad 13-moment approximation thus con- 
tains the implicit assumption that the contributions to E^ terms (k = n) from the higher 
moments are negligible. The validity of this assumption is investigated in the present 
research by a direct comparison of the induced pressure and temperature gradients com- 
puted by the two methods. Since the same collision model and the same power series 
expansion of f e are used in both sets of calculations, the exact results are obviously the 
correct standards in these comparisons. 

A number of additional assumptions are made in the present research to facilitate 
the exact and Everett solutions to the approximate kinetic equations. These assumptions 
are reflected in the following plasma conditions: the heavy particles are infinitely mas- 
sive and at rest relative to the laboratory, the ions are distributed in such a manner as to 
yield zero local charge densities, the applied and induced magnetic fields are zero (the 
latter in strict violation of Maxwell's equations, but consistent with a nonrelativistic 
treatment), all spatial inhomogeneiti.es are induced by the application of a constant elec- 
tric field along the axis of a cylindrical geometry, viscosity is neglected so that the axial 
electron diffusion velocity does not vary radially (through second order), and the system 
is independent of time. 


THE EXACT FIRST-ORDER SOLUTION 


The derivations of plasma transport properties are usually based in some manner 
upon the electron Boltzmann equation (ref. 6) 


at / 1 2kTg J ay 



( 2 ) 


where y is the reduced electron particle velocity defined by 
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(3) 


y = 



and Tg is the spatially homogeneous contribution to the electron temperature. 

If the applied electric field and the electron velocity distribution function are 
expressed through second order as 

E = kE 

and 

/ \ 3//2 

f e = f^ 1 + 4>1 + $2) = n e\T-^5j e " y2 ( 1 + *1 + ^ 2 ) 


(4) 

(5) 


the following simplification of equation (2) is obtained from the application of Meador's 
collision model (ref. 5) and the special assumptions given in the preceding section: 



«(o)W . 


o\l/2 


= " n h r [mf) 1(^1 + ^2 ~ $1 ~ <^ b db de 


( 6 ) 


The only aspect of the collision model which appears explicitly at this stage in 
equation (6) is the assumption that the effective collision partners for electrons can be 
regarded as heavy scattering centers corresponding to a single species of heavy particles 
with number density n^. Corrections for non-Lorentz gases, especially those involving 
electron-electron encounters, are discussed subsequently. 

Since the function <p j is proportional to E by the definition in the preceding 
section, so that the m e cf e moment of f^^i must give the linear form of Ohm's law 
(that is, the proportionality of the axial electron diffusion velocity to E), the assumption 
that the axial electron diffusion velocity does not vary radially through second order pre- 
vents from being x and y dependent in equation (6). The elimination of these 
derivatives and the third-order contributions E ^2 from equation (6) thus yields the 
complete second-order equation 
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( 7 ) 



the first-order form of which is 

, ,1/2 

2eE/ m e 


m e l2kT° 


r-7. - 



* 4>l) b 


db de 


( 8 ) 


Equation (8) can be solved by the assumption that <p^ has the form y z multiplied 
by a function of y. The result is 


_ _ 2eTR 04 E / me \ 1/a (4/D-l 
1 m e R 13 \2kTgy : 

if the Tg-dependent collision time 

mp(r 

r = — — 

e 2ng 


and the collision integral (ref. 5) 



b db de 


r 13 ( m e V^ 2 -4/S 

— ^ 1 n Y Yz 

n h R 0 4 r \2kTgy 


(9) 


( 10 ) 


(ID 


are employed. The Rjj integrals are defined by 



( 12 ) 


As explained in reference 5 and mentioned in the preceding section, the effective 
interaction parameter S in equations (9), (11), and (12) formally appears as the expo- 
nent in a single type of inverse-power electron — heavy-particle interaction potential but 
can be chosen semiempirically to make the corresponding collision model reliable for 
some real plasmas with electron-electron encounters and multiple heavy species included. 
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Only for true Lorentz plasmas (slight ionization or fully ionized gases with large ionic 
charges), however, can £ have physical meaning. 

The first-order electron current density j and the dimensionless residual heat 
flux Q are found from the velocity moments of equation (9) to satisfy 


j = - 


en e Vg = -en 


,2kT^ 

o| e 

m f 


,1/2 


ft) = o-E 


(13) 


and 


Q = ^-|J 0 = .2(|^/^ 

2 0 £m e l 2kT ° 


\l/2 


E - -2(1 + 

* 2^0 “ 


= < 14 > 


if the general nondimensional flux ft is defined by 

ft = <y 21 y) = y2l ^ f e d ^e (15) 

The use of the spatially homogeneous contribution to the electron number density in 
equation (13) is dictated by the restriction to first order at this stage. 


THE EXACT SECOND-ORDER SOLUTION 


Substitution into equation (7) of equation (9) and its velocity derivative 


gives 


hi __ .^f^/M 1/2 r (4/ ?) -3r {y2 _ (£ . 4h | 
dy z m e R 13^ \2kTg / L 


o\V* t 


(2kT 

I- 


-W + y y 9 -^j * m e kT°Cl Vl4/£) ' 3 [ iy2 ‘ + 5 - 4 ) r i 


/2kT, 


, 1/2 


= - n hn 


m f 


§(<t>2 ~ 02l) b 


db de 


(16) 


(17) 


the solution of which provides the exact (through second order) expressions for the 
induced pressure and temperature gradients. Equation (17) differs from the usual 
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second-order formulation of Chapman and Enskog in the explicit appearance of (p 2 on 
the left side. 


By using the collision integral 


I 


„ ,b db de = _ . 

2n hR04 T UkT? 


^13 ( m e \ 1/2 y -4/? 


(3r| - y 2 ) 


(18) 


derived in reference 2 and the definition 


^ (l - cos 2 x)b db 
y ^1 - cos db 


(19) 


one easily confirms that the correct solution of equation (17) is 



where r 2 is x 2 + y 2 . 

In a manner analogous to the development in equation (14), the complete second- 
order formulation of the dimensionless residual heat flux is found from equations (5), 
(14), (15), and (20) to be 


m £ 


Q = h~ 2 Po = 'gpoUro 


1/2 


E • j + j 


m. 


JLi 

2p°[2kT°, 


\ 1/2 


E • j - k- 


U 


@0 


( 21 ) 


Moreover, since the transverse components of the dimensionless electron diffusion veloc- 
ity /3 q are zero from equations (5), (9), (15), and (20), the ordinary heat flux q is 
obtained from equation (21) and the equation of state 


P° = n°kT° 


as 



T + j|E 


5 - 7 + k - - „ + 1} p°v e 


I 
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( 22 ) 



This expression automatically satisfies the energy balance relation 


V • q = E • j 


( 23 ) 


The induced number -density and pressure derivatives are likewise derived from 
equation (20) by the following integrations: 


and 


9ne 

Hx 



d _ _ (£ + 4)n°x / eE 

6 2 £ kT" 


(24) 


!£e = JeL A (VvV*, . . (s - 4 ^/4 

9x 37 j 3/2 9x J 6| [kT°y 


(25) 


whereupon 

2({ +2)Tgx / eE \2 
8 * " n°k 8x ‘ n° 8x ' ' 3 « IfTgJ 

The spatially homogeneous contributions to n e and T e are used in the coefficients of 
the derivatives in equation (26) to avoid terms containing E^. 

Similar expressions (with x replaced by y) hold for the y derivatives, so that 
equations (25) and (26) comprise the exact second-order induced pressure and tempera- 
ture gradients mentioned in the Introduction as a purpose of the present research. 


GRAD 13-MOMENT APPROXIMATION 


The final effort of the present research is the determination of the induced pressure 
and temperature gradients by using the Grad 13-moment approximation and Meador's col- 
lision model (ref. 5) to close out the macroscopic equations of change. Numerical com- 
parisons with equations (25) and (26) should yield valuable insight into the applicability of 
this technique to nonequilibrium plasmas because a completely common framework (i.e., 
basic assumptions and expansions in powers of E) is provided for both methods. 


If the reduced electron particle velocity relative to the electron frame of reference 
is defined by 
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Everett's formulation of the Grad 13-moment velocity distribution function through 
second-order terms can be written as 


f e,Grad - n e(^r] e " u2 j 1 - y( 2u2 “ 3 ) + f f 2 ’ f (?1 " f " 


SL 

p e + — ~(u - 3kk. 


: uu 


( 28 ) 


/3 2 

for the present problem (see the appendix). The term - -^^2u 2 - 3^ within the braces 

in equation (28) and also the form of the tensor product follow from the conversions of the 
electron temperature and the traceless pressure tensor from the electron frame of ref- 
erence to the laboratory frame of reference, 
o 

- 

With P e anticipated to be spatially homogeneous through the second-order level, 
which is subsequently confirmed in equation (35), the pertinent macroscopic equations of 
change are obtained as follows from the application of equation (28) to the m e c e x and 
m e CgC e ^x moments of equation (2): 

= m e J c e,x(^) d ?e = " n h m e J c e (c ejX - c^f e b db de dc e 


4R 13 n e m e^l,x[ 2kT e ] f _ v 2 2 1— (4/ ^ )/ 2 5\ j- 

= ' 5,3/2 R04 ; (=rj T-f 4 ” 

1 /p 

(I - «)Rl3B-l,5ngm c ^ 1 ^2kl"'\ 


5{R§ 4 ^ 


\ m e 


(29) 


and 


a / 2 2 

^In e m e ^c e c| 5 


,x)) = m e J 4 c e,x\^) c dc e = -%m e § c 3 ^c ejX - c^fgb db de dc € 


4R^n°mJi v /2kT 


o\3/2 


L 13 11 e 1 “e H l^ 


J e ^ 2 >|r 3 - (4/l) (y 2 - 1] df 


57t 3 / 2 R 04 t 

= (3£ - 2)(3| - 4)R 13 R. 1 ?5 nQm e ^ 1>x /2kTg\ 3/2 


5| 2 R 2 4 t 


m e 


(30) 
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The successive simplifications of the collision integrals in these expressions cor- 
respond to the use of Meador's collision model (ref. 5), the concepts of inverse collisions 
and microscopic reversibility (ref. 6), and equations (11) and (28). In addition, u is 
replaced with y in the application of f e Grad to the collision integrals, the distinction 
being unimportant here because it represents third-order contributions in view of the 
factor 

The Grad value of /3j x is compelled by the solution of equation (23) to satisfy 

x / m e \ ^ sr . -r _ ^e T Gra.6? i f m e \ ^ ( eE \ 2 /o^ 

2p°\2kT ° ) 2m e \2kTgJ (kTgJ 

the second equality of which follows from equations (10) and (13) with the use of 



R§ 4 (l3£ 2 - 16| + 16) _ 
~«13 R -1,5< 3 « - 2 >* T 


(32) 


from reference 3 to preserve a consistent treatment of the 13-moment approximation. 
Accordingly, equations (29) and (30), respectively, become 



U - 4)(l3| 2 - m + 16)pgx/ eE \ 2 
405 2 (3£ - 2) \kT®y 


(33) 


and 



(3| - 4)(l3| 2 - 165 + 16)nge 2 E 2 x 
205 3 m e 


(34) 


A most important feature of equation (34) is that \c|c| >x ^ is not one of the 
13 moments; consequently, a distribution function must be employed to express this quan- 
tity in terms of the more familiar variables. The error thus introduced by the Grad 
13-moment approximation is expected to be quite large by analogy with the results for 
higher moments calculated in reference 2, and it may even overshadow the errors origi- 
nating from the use of the Grad approximation in the collision integrals. Only the latter 
discrepancies appear in studies of the 13 moments in spatially homogeneous plasmas or 
in the investigation of the induced pressure gradient in the present research. 

Additional emphasis on the error inherent in the closing out of (44,x) with the 
Grad 13 -moment approximation is provided by a study of equations (9) and (20) in the 
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special case of Maxwellian particles (| =4). In particular, since the traceless electron 
pressure tensor is calculated from equation (20) to be 


Pe - n e m e 


f( C e c e> "l( C e) U ) “ "(* 


i i + j j - 2kkJ 


a 4(3{ + 4)R 04 R 22 p°^ 


15A^Rj 3 


(35) 


the substitution of this expression and equations (13), (21), (24), and (26) into the exact 
second-order solution of equations (5), (9), and (20) yields 

fe,exact(l ■= 4) - *e,Grad - 4°VoV 2 - f)( 2 ^ - r 2 ) - - 20y2 + 


^) + ^ 4r4 ' 12y2 + 3 ) 


(36) 


The function f e? Grad is gi yen in equation (28). 

Although the terms in equation (36) which are additional to the Grad 13-moment 
function contribute neither to the direct calculation of the 13 moments nor to the closing 
out of the collision integrals in the corresponding macroscopic equations of change, the 
third term in the bracket does contribute to the closing out of the quantity 
8 (n e m e ^cf c 2 ,x))/ 8x in equation (34). More specifically, the use of equation (36) on the 
left side of equation (34) yields the following second-order result with the aid of equa- 
tions (10) and (13): 


gH n e m e 


c 2 c 2 ' 
^e^e,^ 


8 

f 2 

8x 

m ej c e ( 

8 

5kp e T e 

8x 

m e 


27r 3 / 2 m f 


+ P 3 / 2 t 2 I e ~ y2y x( 4 ^ 4 “ 12y2 + 3 ) d r 

^'e“'> / 2 y^ 4 y 4 _ 12 y 2 + 3) dy 


5kTe 

m e 


n e k 


8Te 

9x 


^Pe 
+ 8x 




5kTg| 
m f 


0l . 9T e 


npk 


9 P e 


+ • 


8x 8x 


^Pe x / eE 
5 VkT° 


(37) 
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Since the pressure gradient is zero from equation (33) when £ = 4, the combination 
of equation (37) with equation (34) yields the following complete £ = 4 temperature 
derivative: 


8T e 

8 x 




(38) 


The first and second numerical fractions on the right side of this expression correspond, 
respectively, to the contribution from the collision integral (which is the same for both 
the Grad and the exact distribution functions if £ = 4) and the contribution to 
8^n e m e ^Cgclpj^/sx from the last term in equation (36). In addition, the total value for 

STg^Sx agrees with the exact result in equation (26) when £ = 4; hence, this analysis 
firmly establishes the importance of one of the terms neglected by the Grad 13 -moment 
approximation. 

Another interpretation of the last term in equation (36) is obtained from the following 
heat flux expression deduced from equations (30) and (37): 


^x 



8Tc 


C l-Sf + C 2 


1ST + °3 e2x 


(39) 


Except for the coefficient C3, which is zero in the 13-moment approximation but is quite 
significant in the exact solution, equation (39) is the familiar transport relation with 
transport coefficients Cj and C2. Since the 13 -moment values of q x , Cj, and C2 
are fairly accurate, large errors in the temperature or pressure gradients (or both) must 
exist in that approximation to compensate for the serious error in C3. 


Finally, the use of only the Grad 13 -moment distribution function to perform the 
indicated average on the left side of equation (34) gives 


_o, 8T e 9 Pe ( 3 ? - 4 )(l 3| 2 - 16 £ + 16 )pgx/ eE \ 2 
ne 8 x + 8 x " 10 0£3 \kTo/ 

so that 


(l 3| 2 - 16 £ + 16 ) 

I 2 T°x 

( eE \ 2 

200 £ 3 ( 3 £ - 2 ) ! 

fTg) 


( 40 ) 


(41) 


with the aid of equation (33). Equations (33) and (41) thus comprise the Grad 13-moment 
values of the induced pressure and temperature gradients mentioned in the Introduction 
as a purpose of the present research. 
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COMPARISON OF RESULTS 


A convenient representation of the preceding results is through the following defi- 


nitions of nondimensional pressure and temperature derivatives: 

X = 1 A T ef d Pe 

rp°\ eE / dr 


and 


Y = - 


rT 


fkT°\ dT e 
,o\eE / dr 


(42) 


(43) 


Accordingly, from equations (25), (26), (33), and (41), 


X exact 


= LiA 

6 £ 


x Grad 


(| - 4)(l3| 2 - 16| + 16) 


40r(3^ - 2) 


(44) 

(45) 


Yexact^^^f^ 


and 


Y Grad = 


_ (l3g 2 - 16$ + 16) 2 


200£ 3 (3£ - 2) 


(46) 


(47) 


Numerical calculations appropriate to equations (44) to (47) are presented in table I 
for a variety of effective interparticle interaction potentials ranging from the fully ionized 
Lorentz plasma (ij = 1) to a gas of rigid spheres (£ = °°). The expectations discussed in 
the paragraph following equation (34) are clearly realized: Although the 13-moment val- 
ues for the induced temperature gradient are quite poor over the entire range of £ 
(Maxwellian particles included), the errors in the induced pressure gradient are far more 
tolerable except for very soft or very hard force laws. 

A PROPOSED USE OF THE EXACT SOLUTION 

o 

The comparisons given in the present research are between the exact (through E 
terms) and 13 -moment solutions of an approximate Boltzmann equation. Although the 
collision model employed is exact for Lorentz plasmas and has been shown to be valid 
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(ref. 5) in the first-order theory of real fully ionized gases (electron-electron collisions 
included), its applicability to second- and higher-order treatments of non-Lorentz plas- 
mas has not been demonstrated. Consequently, the present exact results should not be 
regarded as the final answer in the non-Lorentz case but rather as an indicator of the 
magnitudes and sources of errors inherent in the 13 -moment approximation. 

More specifically, the importance in the exact solution (given by eq. (36)) of the 
terms which are additional to the 13-moment distribution function requires that a mini- 
mum of four new moments must be added to the original 13 moments to complete the 
moment representation of equation (36) and thereby derive a minimally acceptable Grad- 
like approximation for the present problem. This approximation can then be employed 
in the usual Everett manner to close out the macroscopic equations of change corre- 
sponding to the rigorous collision integrals instead of the present collision model. Since 
the exact solution of the approximate Boltzmann equation serves in this method only to 
establish the minimum number and types of moments to be incorporated in the Grad-like 
approximation, the limitations on the original collision model should have a trivial impact 
on the final results of the ensuing closing-out procedure. 

The additional terms in equation (36) can, of course, be represented by more than 
four moments, in which case the proposed closing-out procedure will give improved 
results because the increased number of macroscopic equations of change will offer the 
rigorous collision integrals a greater opportunity to exert an influence. 

Another important advantage of this proposed use of exact solutions concerns the 
forms of the additional terms in equation (36). If moments were added in the usual Grad 
fashion (employing Hermite polynomials) to improve the 13 -moment approximation for 
the present problem, the convergence would be extremely slow because of the difficulty 
in describing by this method the y~% dependence of one of the terms. The proposed 
technique introduces the y ~ 2 dependence directly, so that a very rapid convergence is 
expected. Further improvements in the proposed Grad-like approximation can be 
obtained by using a more realistic | in the exact solution, as opposed to equation (36) 
which is based on Maxwellian particles. This alteration may, in fact, be necessary if 
some of the important y dependencies in the exact solution are functions of the inter- 
action potential. 


CONCLUDING REMARKS 

Calculations through second order in the electric field have indicated that large 
errors can occur when the Grad 13 -moment velocity distribution function is used to close 
out the macroscopic equations of change for a steady-state plasma. Much of the difficulty 
in computing the induced temperature gradient, the values of which are especially poor. 
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can be traced to the failure of the 13-moment approximation to predict accurately an 
important higher moment appearing in the macroscopic equation for the heat flux. The 
pressure gradient, on the other hand, does not involve the closing out of such higher 
moments, so that the 13 -moment predictions of this property are in much better agree- 
ment with the exact second-order values. 

The present results, and those of previous research, suggest that some of the past 
studies employing the Grad 13-moment approximation perhaps should be reevaluated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 19, 1970. 
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APPENDIX 


TRANSFORMATION OF THE GRAD 13-MOMENT VELOCITY DISTRIBUTION FUNCTION 


The purpose of this appendix is to derive equation (28) from the original Everett 
formulation (ref. 4) of the Grad 13 -moment approximation, which is written as 


^e,Grad - n e 


m. 


1 3/2 


27rkTi 


expl 


m e 

2kT 


— \2 


r(c e - v e ) 



m e u. -* \2 


Cp - V, 


2kTp V e ,e 



1 m e 

Pe 2kT^ 


o 





Note that the primed quantities are measured relative to the electron frame of reference 
and that n e is the same in both the electron and the laboratory coordinate systems. 

The following relations are fundamental to the transformation from equation (Al) to 
equation (28) if terms higher than second order are neglected and if the differences 
n e - ng and T e - T° are proportional to E 2 for v|] as in the present problem: 




and 


Pe = n e m e 


((je ~ Ve)(c e - v e )) - |^(c e - v e ) 2 )lf = n e m e ((c e c e y 


v e v e 



+ (Pe - P4)U - 2pg^kk * P e 



Equation (A5) anticipates the radial electron diffusion velocity to be zero. 


p^U 


(A5) 
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APPENDIX - Concluded 


o_ 

If one further anticipates that P e is second order, which is confirmed by the 
exact result of equation (35), equations (Al), (A4), and (A5) combine to yield 



Since T e differs from by a second-order term, the introduction of iT of equa- 

tion (27) in the non-Maxwellian terms of equation (A6) is consistent with a second-order 
analysis because the vector and tensor factors are first order or higher. 

Equation (28) is immediately obtained from equation (A6) when equation (A2) is 
employed to rewrite the Maxwellian contribution in the form 
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TABLE I.- NONDIMENSIONAL INDUCED PRESSURE AND TEMPERATURE GRADIENTS 



X 


Y 


13 moment 

Exact 

13 moment 

Exact 

1 

-0.975 

-0.500 

0.845 

2.000 

2 

-.113 

-.167 

.203 

1.333 

4 

.000 

.000 

.200 

1.000 

CO 

.108 

.167 

.282 

.667 


NASA-I.angley, 1970 25 L-6877 
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